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A viscous compressible model of soap film flow and its equivalence with the 

Navier- Stokes equations 

Petri Fast* 

Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore, CA 

(Dated: October 24, 2005) 

We present a quasi-two dimensional model of flowing soap films that bears striking similarity to 
the compressible Navier-Stokes equations. The variation in soap film thickness that is commonly 
used for flow visualization in experiments is analogous to density variations in the Navier-Stokes 

equations. When the soap film flow velocity is comparable to the Marangoni elastic wave velocity we 
recover the compressible Navier-Stokes equations and the soap film behaves like a two-dimensional 
isothermal viscous gas. 
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Fast flowing soap films have been used extensively 
as an experimental realization of two-dimensional fluid 
dynamics [1]. The purpose of this Letter is to present a 
new quasi-two dimensional model of flowing soap films 
that bears striking similarity to the compressible Navier- 
Stokes equations. In a particular case, we can for- 
mally identiiy the viscous soap film model with the two- 
dimensional compressible Navier-Stokes equations: The 
variation in soap film thickness that is commonly used for 
flow visualization in experiments is analogous to density 
variations in the Navier-Stokes equations in a quasi-two 
dimensional setting. The effective viscosity in our model 
is variable and depends on the film thickness. 

The original papers [2, 3] lay out the experimental and 
theoretical foundation for considering flow in flat soap 
films as classical two-dimensional flow. A new experi- 
mental technique developed in the late 1990's[l] further 
popularized this approach by enabling very large, long- 
time stable, gravity driven soap films that have been used 
for studies of two-dimensional turbulence [4-6], fiuid- 
structure interactions [7] and Shockwave dynamics [8]. 

Most previous theoretical work on modeling fiowing 
variable thickness soap films has considered them as rep- 
resented by the two-dimensional constant density incom- 
pressible Navier-Stokes equations[2]. However, Vorobi- 
eff et aL[4] have pointed out in a series of papers sig- 
nificant thickness variations in rapid soap film flows. 
Chomaz [9] derived recently a new model of soap film 
flow and argued that: (1) The model in an inviscid limit 
is equivalent to the two-dimensional compressible Euler 
equations, and (2) no correspondance with the compress- 
ible Navier-Stokes equations could be found. 

In this Letter, we establish a connection between the 
compressible Navier-Stokes equations and a model of 
flowing soap films. We derive a quasi-two dimensional 
viscous compressible model of fast flowing soap film 
that is applicable in subsonic and supersonic regimes. 
The model is derived systematically from the three- 
dimensional incompressible Navier-Stokes equations cou- 
pled with a model of surfactant transport and capillary 



forces on the free surfaces. 

We consider a flowing soap film of thickness H and 
width L and scale z ^ H and the lateral directions a;, y ~ 
L. In typical experiments H « lO/itm and L « 10cm. We 
develop a thin film model using the small parameter e = 
H/ L. The film is driven by an incoming flow rate Q which 
determines a characteristic flow velocity U — Q/HL. An 
important characteristic of soap fllms is the Marangoni 
stress that arises from variations of the surface tension a 
that depend on the local surfactant concentration V on 
the free surface z = h{x, y, t). We assume a linear surface 
tension model[10] cr = CTo — fTrF with parameters Ua, (Jr- 
We scale the velocity u,v ^ U, w eU, and the pressure 
p ~ ecJrn/L to balance the capillary forces, and define 
the mean surface tension <jm= (^a — (^r^m, and the mean 
surfactant concentration r^. 

The scaling leads to five nondimensional groups: the 
Reynolds number Re = pUL/fiQ, the elastic (Marangoni) 
Mach number Ma = U{pH/{arTjn)y^^ , the bending 
(capillary) Mach number Mb = e~'^U {pH / (7^)^/"^ , the 
surface Schmidt number Sg = l^o/{pDsurf ) for the diffu- 
sion of r, and the bulk Schmidt number Sc = fio/ipDc) 
for the diffusion of c. We define relaxation constants 
A, K* that describe the exchange of surfactant concentra- 
tion r at the surface with the surfactant concentration c 
in the interstitial fluid. 

A new formulation of two-dimensional viscous com- 
pressible soap film flow is given by 
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based on the primitive variable equations introduced by 
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Chomaz[9]. Here / is the identity tensor. A full deriva- 
tion of Eqs. (l)-(4) is provided below after we discuss 
some of the properties of this model. 

A key result of this Letter is the new formulation of the 
viscous stress tensor //(Vu + Vu"^) + 2r]"{V ■ u)I in the 
soap film model: This is formally identical to the viscous 
terms in the compressible Navicr- Stokes equations with 
dynamic viscosity rj = 2/i/Re and dilatational viscosity 
T]" = h/Bje. Actually, the viscous terms in Chomaz [9, 
Eqs. (3.14)], andlda&Miksis [11, Eqs. (14) (18)], which 
were written out in nonconservative form, are identical to 
the viscous term in Eq. (1). The benefit of rewriting the 
viscous terms in conservative form is that it reveals the 
similarity with the compressible Navier-Stokes equations. 
Past work[9, Sec. 4.2 on p. 404] expressed doubts that a 
soap film model can be compared with the compressible 
Navier-Stokes equations since it appeared that the soap 
film viscous stress could not be expressed in a form such 
as Eq. (1). 

In fact, Eqs. (l)-(4) are formally equivalent to the com- 
pressible Navier-Stokes equations in flows where disper- 
sive effects are insignificant (1/Mb 0) and surfactants 
are insoluble (A 0). The thickness h and concentra- 
tion r in the soap film model correspond to the density 
and pressure, respectively, in the Navier-Stokes equations 
with a special equation-of-state. We consider two special 
cases. 

First, we consider a finite Reynolds number version 
of the "inviscid supersonic soap film" limit in Ref. [9, 
Sec. 4.2]. The interstitial soap concentration c = Cq is 
assumed to be constant. Eqs. (l)-(4) simplify to 
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where the dynamic and dilatational viscosity is rj = h. 
The ratio T/h is conserved along streamlines in this 
limit [9] which corresponds to the behavior of p/p in 
isothermal flows. Therefore, the soap fllm Eqs. (5)-(7) 
are identical to the compressible Navier-Stokes equations 
for isothermal fluids without viscous heating. In the in- 
viscid limit Re — > oo these equations are identical to 
Eq. (4.13) in Ref. [9]. 

Second, we consider incompressible flow in a variable 
thickness soap film. Eqs. (5)-(7) simplify to 

^ + V • {huu) = -V<A + • {,?(Vu + Vu^)} , 

ht + V- (hu) = 0, V • u = 0, where we have define a scaled 

surfactant concentration (f) = T/Ma^. This model corre- 
sponds formally to the incompressible variable density 



Navier-Stokes equations with a density dependent vis- 
cosity, and with acting as a pressure-like variable. The 
inviscid limit Re ^ oo reduces to Ref. [9, Sec. 4.1.2.]. 

We note that air drag effects can significantly change 
soap film flow measurements over flow distances of tens of 
centimeters [12, 13]. However, recent experiments operate 
in a regime where air drag effects are secondary and the 
dominant forces are due to the dynamics of the thin liquid 
layer[14]. Future work will consider air drag effects[12, 
13] in the model (l)-(4). 

We discuss an "acoustic" limit as the simplest flow that 
can be used to illustrate compressibility effects. Consider 
the small perturbations u = + m, h ^ 1 + h, and 
r = 1 -|- r on top of a quiescent base flow u, h, and f . 
The leading order linear equation is 
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The last term in Eq. (8) causes viscous dissipation of 
the wave solutions and is omitted for simplicity in the 
following discussion. The dispersion relation for solutions 
of the form u = exp(i(wi + kx)) reveals two types of 
limiting behavior. Pure bending waves (for Ma, Re oo) 
are described by 



1 



Mb^ dx* 



and the dispersion relation lu = ±fc^/Mb. These waves 
are dispersive [15] with a group velocity cLu/dk — 2A;/Mb 
that is twice the phase velocity. This is similar to the 
dispersive effect of the capillary terms in water wave 
theory [15]. The Navier-Stokes or Euler equations do not 
exhibit this type of behavior. 

Marangoni waves arise from the limit Mb, Re — > oo, 
where the small amplitude dynamics is governed by 



1 d'^u 



These elastic waves arise from stretching of the film 
that induces tangential stresses to restore the equilib- 
rium concentration. Marangoni waves are analogous to 
sound waves in gas dynamics with the dispersion relation 
LO = ±fc/Ma, and a (nondimensional) phase and group 
velocity 1/Ma that is independent of wavenumber. Slow 
flowing soap films can be viewed as 2-d incompressible, 
when the flow speed is a small fraction of the elastic sound 
speed (Ma <C 1). However, many experiments operate in 
a regime where Ma = 0(1) and 2-d compressibility ef- 
fects are clearly visible[16]. Marangoni waves have been 
measured in recent soap film experiments [8, Table 1]. 

We present the details of the derivation for complete- 
ness and to simplify some aspects of past work [9, 11]. 
The dynamics of the soap film is governed by the three- 
dimension incompressible Navier-Stokes equations with 
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constant density p and viscosity /i. Wc define two- 
dimensional quantities u = {u,v), Vp = {pxiPy)i and 
V^M = Uxx+Uyy and use subscripts to denote derivatives. 
Wc assume the soap film flow is symmetric about the cen- 
ter surface and impose Uz = 0, = 0,^ = 0, = 
at 2; = 0. Hence we ignore the bending mode and are 
only considering the symmetric (peristaltic) mode [2, 9]. 
The thin layer of fluid is surrounded by free surfaces 
z = ±h{x, y, t) on which surfactant is transported. The 
surfactant concentration c{x,y,z,t) in the interstitial 
fluid satisfies [9, 17] an advection-diffusion equation with 
a diffusivity Dc and a free-surface boundary condition 
Dcdc/dn = j. The flux j = {Kc — T)/t models surfac- 
tant transport between the surface and the interior of the 
film with transport parameters K, r. The use of both an 
interior and a surface concentrations c and F in a soap 
film model is due to Chomaz[9]. 

The surfactant on the on the free surface z = h{x,y,t) 
has a variable concentration T that is transported accord- 
ing to [18] 

BT 

— + V, • (Fu) + {2M){n ■ u)r = V^F +j (9) 

where the flow and derivative operators are restricted 

to the surface, the flux j allows exchange of surfactant 
with the interior concentration c, and 2A4 is the mean 
curvature (see [19, 20] for details). We denote the average 
concentration Cm, which yields the average surfactant 
concentration F^ = KCm- We scale A = L/{tU), and 
K* = K/H. 

The nondimensional bulk equations are, on dropping 

the primes, 
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and the divergence free condition V • u + dw/dz = 0. 

The free surface z = h{x,y,t) satisfies a kinematic 
boundary condition ht + uhx + vhy = w. The vis- 
cous surface stresses are balanced by capillary forces and 
Marangoni stresses arising from the surfactant concen- 
tration distribution. The nondimensional surface tension 
is given by cr = 1 + C (1 — F) , where C = f rFm/cm = 
e^Re/Ma^. The scaling C = ©(e^) is determined by a 
distinguished limit arising in the tangential stress con- 
dition (see below). The normal stress condition \—p + 
2/xon • D • n] = 2gM. in nondimensional form is 



{p + 2M (i + c(i-r))) 

-e^{hlux + hlvy 



2MV 



{hxUz + hyVz + V • u 



Re^^ 

+ hxhy{Uy+Vx) 

- hxWx - hyWy)} + 0(e''), 



where = l + e^\SIh\'^ and 2M = V^h + 0{£^)[lQ, 20]. 

The tangential stress conditions [s^ • D • n] = ■ VgCr 
and • D • n] = • VgCr in dimensionless form are 

Uz + £^[-2hxUx - hy{uy + Vx) +Wx- hluz 



hxhyVz + 2hli 



Vz + e'^[-2hyVy - hx{Uy -I- Vx) +Wy - kyV z 

Re 

- hxhyUz + 2hlwz] = ~^^^^^v + Oi^'^)- 

Note that we have the choice of at least two scalings of the 
coeflacient C- In the first case, C = e^Mb^/Ma^ = ©(e^), 
and the surfactant concentration dependent terms do 
not enter the dynamics through the normal stress con- 
dition (10) at leading order which leads to the soap film 
model (l)-(4). This is the distinguished limit considered 
in Chomaz[9], following earlier work [10, 11]. In the sec- 
ond case, C = 0(1), the leading order fluid velocity would 
be slaved to the surfactant concentration[10]. This is not 
a physically relevant scaling for gravity driven fast soap 
fllm flows. The C = 0(1) scaling is not considered here 
any further. 

We expand the solution q = {h,u,v, F, c, p) = q-l-e^q-i- 
0(£^) in £. We will obtain equations of motion depending 
only on {x, y, t) by averaging across the thin layer in the 
z direction. 

The leading order tangential stress condition = 
implies u = u(a;,y,t). The vertical velocity w satisfies 
w = — zV • u. The leading order velocity ti = \i{x, y, t) is 
determined by the 0{£^) equations as in the derivation of 
thin jet models[10, 21]. The normal stress condition (10) 
at leading order is 

^oj 2Mb2 ^ 

The leading order mass conservation and surfactant 
transport equations are Cz = 0, 



ht + uhx + vhy = w = — /iV • u, 
ft + ufx+vfy + f{V-vi) = A(c-f) 



SgRe 



V^F. 



The interstitial concentration satisfies Cz = 

— \e^RcScK*{c — F), which allows integration of 
the leading order to obtain Eq. (4). The momentum 
equations at 0{e^) are 

and the tangential stress conditions are 
Re 



Uz = 



Vz = - 



^^2^x + 2hxUx + hy{Uy + Vx) 

+ h{V ■ u)x + 2hxV ■ u. 

Re - - _ - _ 

■— ^Fj^ + 2hyVy + hx{Uy + Vx) 

Ma 

h h{V ■ u)y + 2hyV ■ u. 
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The leading order momentum equation is rewritten in 

divergence form as 



Ref;.^ + 
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= dx {hvx + huy) + 2dy {hvy) + 2dy (/iV • u) 

Eq. (1) follows since we can define the rate-of-strain ten- 
sor D = (Vu + Vu'^)/2 and use 2V • (/iD) to simphfy 
the RHS. This completes the derivation. 

In summary, we have presented a quasi-two dimen- 
sional viscous compressible model of soap film flow that 
shares many similarities with the Navier-Stokes equations 
but that also presents some intriguing differences. Future 
research is needed to determine the appropriate bound- 
ary conditions for the model since the analysis herein and 
in past work [9] is only valid away from physical bound- 
aries extending through the film. Further, the vorticity 
w = Vx — Uy in a soap film satisfies, ignoring molecular 
diffusion (Re — * oo), 

|:+u.V. = -.(V.u) + ^ 



v/i X vr. (10) 



The baroclinic[22] term Vh x VP 

as a vorticity source in regions where the isocontours of 
thickness and surfactant concentration do not coincide. 
This effect is absent from constant density Navier-Stokes 
flow. Numerical simulations using Eqs. (1) (4) will be 
necessary to investigate the importance of baroclinic ef- 
fects, surface bending elasticity contributions, viscosity 
and compressibility in fast flowing soap film experiments. 
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University of California Lawrence Livermore National 
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SUPPLEMENT: 

A viscous compressible model of soap film flow and its 
equivalence with the Navier- Stokes equations 



SOAP FILM SURFACE GEOMETRY 



In this section, wc summarize differential geometry results that arc needed for describing the surface geometry of 
the soap film. We derive expressions for the mean curvature of the surface following Struik[20, Sec. 2-7]. We consider 
the dynamics of a thin fluid layer that is symmetric about the center surface 2; = of the film. 

The surface of the soap film is located at z = ±h{x, y, t). In the following, the time dependence is implied and left 
out from the notation. All the expressions are dimensional. (Note that in the non-dimensional, scaled, setting the 
surface height gradients become V/i eV/i.) 

We denote the surface vector by ^(a;, y) = {x, y, h{x, y)). The surface normal is 



Ca; X _ ( ^x, fly, 1) 



where s = {1 + + hyY^^. The mean curvature is, following Struik's notation[20], 

EG- 

1 V^/i -I- hlhyy - 2hxhyhxy + h^h^x 



and 



EG-F' = {l + hl){l + hl)-hlhl = s', 



Eg-2fF + eG = {{1 + hl)hyy - 2hxhyKy + {I + hl)Kx) / s. 
For an almost flat surface we obtain 2A^ ~ V^h. 



(11) 



Si - 777172' ^2 - y^^l 



(l + /,2+/j2)l/2' 

where = (1, 0, hx), = (0, 1, hy), and l^j, x ^j^|2 = \-\-h^-\- hy. The surface tangent vectors are 

(1,0, M g (0, 1, f ly) 

The mean curvature is determined from the first and second fundamental form [20] 

/ = \^x\^dx^ + 2^,-^ydxdy+\^y\''dy^, (13) 
II = ^xx dx"^ + 2ixy ■ n dx dy + S,yy -ndy"^, (14) 

where j^xP = 1 + l^yP = ^ ~^ ■ £,y — fix fiy The fundamental forms simplify to 

I = {I + hi) dx'^ + 2tixfiy dx dx + {1 + til) dy'^ , 

II = ^dx-' + '^dxdy + ^dy-', 
s s s 



1 Eg ^ 2fF + eG 

^ = 2 EG-F^ (''^ 

-Oh h h -L h'^h V 

(16) 



2 (1 + |V/i|2)'/2 

Here 

E=l + fil e=- - (17) 

{l + hl + filf^ 

F = flxfly, /=- (18) 

G=l + fil 9=- ^^^-171- (19) 

{1 + hl + filf' 
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VISCOUS SURFACE STRESSES 

We summarize here the dimensional viscous surface stresses at 0(1) and 0(e^)- Define = 1 + |eV/ip. First, recall 
the tangent vectors = (1, 0, ehj.)/^, = (0, 1, ehy)/^, and the normal vector n = x = {—ehx, shy, 1)/^. The 
viscous surface stresses are, to 0(£'*), 

n ■ D • n = {-hxUz - hyV^ - V ■ u + {hlu^ + h^Vy + hxhy{uy + v^) - Kw^ - hyWy]) , (20) 

• D • n = {uz + {-2hxUx - hy{uy + Vx) + Wx - hluz - hxhyVz + 2hlwz}) , (21) 

• D • n = 2^j^^2 + {-'^K'^y - hx{uy + Vx) +Wy - h^Vz - hxhyU^ + 2h'^yWz]) . (22) 



